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ABSTRACT 

We present high resolution X-ray spectra of 14 putative cooling-flow clusters of galaxies obtained with the 
Reflection Grating Spectrometer on XMM-Newton. The clusters in the sample span a large range of temperatures 
and mass deposition rates. Various of these spectra exhibit line emission from O VIII, Ne X, Mg XII & XI, Al 
XIII & XII, Si XIV & XIII, N VII, and C VI as well as all Fe L ions. The spectra exhibit strong emission from 
cool plasma at just below the ambient temperature, 7b, down to 7b/2, but also exhibit a severe deficit of emission, 
relative to the predictions of the isobaric cooling-flow model at lower temperatures (< 7b/3). In addition, the best- 
resolved spectra show emission throughout the entire X-ray temperature range, but increasingly less emission at 
lower temperatures than the cooling-flow model would predict. 

These results are difficult to reconcile with simple prescriptions for distorting the emission measure distribu- 
tion, e.g. by including additional heating or rapid cooling terms. We enumerate some theoretical difficulties in 
understanding the soft X-ray spectra of cooling-flows independent of the classic problem of the failure to detect 
the cooling-flow sink. Empirically, the differential luminosity distribution is consistent with being proportional 
to the temperature to the power of « 1 to 2, instead of being independent of the temperature, as expected in the 
standard multi-phase model. The primary differences in the observed low temperature spectra are ascribed to 
differences in the ambient temperature. 

Subject headings: Clusters: general- Intergalactic medium-Galaxies: abundances-Galaxies: cooling 
flows-Methods: data analysis- X-rays: galaxies: clusters 



1. INTRODUCTION 

It was recognized many years ago that the cores of clusters of 
galaxies have sufficient X-ray luminosity to cool 10 to 1000 so- 
lar masses of X-ray emitting plasma every year (e.g. Fabian & 
Nulsen 1977, Cowie & Binney 1977, Fabian 1994). The details 
of the cooling process are still debated, however, in most mod- 
els, parcels of cooling plasma collect at the center of the cluster 
(Nulsen 1986), forming what is referred to as a cooling-flow. 

This basic picture has remained controversial, however, and 
a number of observational discrepancies and alternate theoret- 
ical interpretations have been raised in the literature. While it 
is well-established that cores of clusters have short (< 10 9 yr) 
cooling times (e.g. White, Jones, & Forman 1997, Peres et al. 
1998, Allen 2000) and cluster cores have been shown conclu- 
sively to contain much lower X-ray temperature plasma than the 
ambient hot outer regions (Canizares et al. 1979, Canizares et 
al. 1982, Mushotzky & Szymkowiak 1988), there is no consen- 
sus about how much mass has cooled from X-ray temperatures. 
Some evidence exists for copious amounts of cooler gas emit- 
ting in the UV (Oegerle et al. 2001), and in Ha emission (Heck- 
man et al. 1989, Crawford et al. 1999). In addition, molecular 
hydrogen has been detected (Jaffe & Bremer 1997, Donahue et 
al. 2000), and the existence of significant quantities of dust has 
been inferred from infrared emission (Edge et al. 1999, Irwin et 
al. 2001, Allen et al. 2001). However, HI absorption measure- 



ments have found no evidence for cold condensed clouds (e.g. 
O'Dea et al. 1998), and while CO emission has been detected 
(Edge 2001) the amount of cold molecular material is still a fac- 
tor of 10 below what is predicted by X-ray cooling estimates. 
Therefore, the connection between these observations and the 
X-ray data is still unclear, and the exact quantity and location 
of cooling-flow byproducts is unresolved. We refer to this as 
the classic cooling-flow problem. 

The bulk of the thermal energy of the cooling intraclus- 
ter medium is thought to radiate at X-ray wavelengths, and 
thus studying the X-ray spectrum is critical to understanding 
cooling-flows and testing cooling-flow models. If the plasma 
cools homogeneously, the density profile of the cores of clusters 
should be much steeper than observed (Johnstone et al. 1992), 
which has led to the conclusion that the cooling plasma must 
condense locally into smaller clouds distributed over a large 
volume (tens of kpc), i.e. in a multi-phase medium. How- 
ever, even ignoring the details of the resulting spatial distri- 
bution, simple thermodynamic arguments show that the inte- 
grated X-ray spectrum of such a cooling flow can be robustly 
predicted. If the blobs of plasma cool in thermal isolation at 
constant pressure, and the dominant energy loss mechanism is 
via X-radiation, then the differential luminosity distribution, i.e. 
the luminosity radiated per unit temperature interval, must be 
proportional to the mass deposition rate, M: 
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dh _ 5 Mk 

dT 2 [rnip 

where k is Boltzmann's constant, and [im p is the mean mass 
per particle. The only free parameter in this expression is M, 
which can be estimated from the density distribution inferred 
from the X-ray image of the cluster core. The resulting spec- 
trum can be calculated using a collisional equilibrium spectral 
synthesis model with an assumed set of elemental abundances, 
and normalizing the contribution in each temperature interval as 
given in Equation (1). To a good approximation, a cluster spec- 
trum should consist of two components: 1) the cooling-flow 
spectrum, as described above, and 2) an isothermal spectrum 
evaluated at the temperature of the background cluster gas. 

Data acquired by the Reflection Grating Spectrometer (RGS) 
on XMM-Newton enable this robust spectral prediction to be 
quantitatively tested for the first time. Surprisingly, the ob- 
served spectra reveal a remarkable systematic deficit of emis- 
sion at the lowest temperatures, as compared to the multi-phase 
model (Peterson et al. 2001, Tamura et al. 2001a, Kaastra et 
al. 2001, Tamura et al. 2001b, Xu et al. 2002, Sakelliou et 
al. 2002). This result has been confirmed through medium 
resolution spatially-resolved spectroscopy using the XMM- 
Newton European Photon Imaging Cameras (EPIC) and Chan- 
dra ACIS observations, where spectral fits have yielded signifi- 
cantly smaller M's than expected (David et al. 2001, Bohringer 
et al. 2001, Molendi & Pizzolato 2001, Schmidt, Allen, & 
Fabian 2001, Ettori et al. 2002, Johnstone et al. 2002). We re- 
fer to the observed deficit of the predicted soft X-ray emission 
as the soft X-ray cooling-flow problem. As we discuss in this 
paper, the soft X-ray cooling-flow problem may or may not be 
related to the classic cooling-flow problem. 

Here we present a systematic study of 14 clusters with the 
RGS to quantify in detail the soft X-ray cooling-flow prob- 
lem. We present four basic results: 1) We find no evidence 
for X-ray absorption by cold material, which had been inferred 
from lower spectral and angular resolution data (White et al. 
1991). 2) We find a severe lack of emission from lower temper- 
ature ions expected in the standard cooling-flow model. 3) We 
demonstrate the ubiquity of significant plasma just below the 
ambient gas temperature, To, down to To/2 in roughly the pre- 
dicted amount. This last result is the most perplexing and is dif- 
ficult to reconcile with proposed explanations for the cooling- 
flow problems. 4) We offer an empirical parameterization of 
the temperature distribution which is consistent with the entire 
sample. 

The paper is organized as follows. In §2 we discuss the ex- 
pected spectrum of the isobaric radiative multi -phase model and 
various associated temperature diagnostics. In §3 we discuss 
the spectral capabilities of the RGS for sources of moderate ex- 
tent. In §4 we describe our analysis methods and present the 
high resolution spectra. In §5 we discuss several qualitative as- 
pects of the spectra. In §6 we describe the model used to set 
cooling luminosity limits. In §7 we present the results. In §8 
we discuss the results with respect to potential modifications of 
the cooling-flow scenario. 

2. DIAGNOSTICS OF THE ISOBARIC RADIATIVE 
MULTI-PHASE MODEL 

It has been argued that the intracluster medium (ICM) should 
remain in collisional ionization equilibrium, even as it cools. 
Both the recombination time scale and electron equilibriation 



time scale through Coulomb collisions, are much less than 
the cooling time scale, evaluated at characteristic X-ray tem- 
peratures and densities in the ICM (Edgar & Chevalier 1986, 
Hicks & Canizares 2001). At high temperatures, the X-ray 
spectrum is produced through bremsstrahlung and line emis- 
sion from abundant elements. The spectra are dominated by 
lines from hydrogen-like and helium-like ions, and iron L shell 
emission. In principle, the temperature distribution can be de- 
termined by three different methods: 1) through the shape of 
the exponential cut-off of the bremsstrahlung spectrum, 2) the 
ratios of hydrogen-like to helium-like lines, and 3) through the 
distribution of line emission from iron L-shell ions. However, 
the bremsstrahlung spectrum constrains only the highest tem- 
peratures present in the spectrum, and the ratio of hydrogen- 
like to helium-like lines is relatively unconstraining since both 
charge states are present over a very broad range of tempera- 
tures. Therefore, Fe L emission is the most useful diagnostic, 
and provides the primary method for determining the tempera- 
ture distribution presented in this paper. 

The isobaric multi-phase model yields a unique spectral sig- 
nature given an isothermal outer cluster temperature and a set of 
assumed abundances. This is shown in Figure 1 for a maximum 
temperature of 8 keV and for 1/3 solar abundances. In Figure 

2, we divide the temperature distribution below 6 keV into four 
bands. The curves are the spectra produced by emission from 
the multi -phase model in bins of 3 to 6, 1 .5 to 3, 0.75 to 1 .5, and 

0. 375 to 0.75 keV. Important line blends are compiled in Table 

1 . Several important diagnostics and aspects of the model can 
be noted. The first is that line emission from ions between Fe 
XXIV and Fe XVII provide tight constraints on the temperature 
distribution. There is roughly comparable emission in the emis- 
sion line blends from each Fe L ion. Also note that Fe XVII line 
emission is predicted to be three times as strong as Fe XVIII. 
O VIII and other hydrogen-like ions are produced at a range of 
temperatures, and thus do not provide strong constraints on the 
temperature distribution. The overall normalization is propor- 
tional to the mass deposition rate, and the spectrum should be a 
strict superposition of all the individual temperatures up to the 
ambient temperature. 

3. THE USE OF THE REFLECTION GRATING SPECTROMETERS 

FOR OBSERVATIONS OF EXTENDED SOURCES 

In order to test aspects of this model in detail, the Fe L spec- 
trum needs to be resolved for a spatially-extended object. This 
could not be accomplished with previous instrumentation; how- 
ever, the launch of the XMM-Newton RGS experiment now 
makes it possible. For a full description of the RGS experment 
see den Herder et al. (2001). The two RGS spectrometers have 
160 cm 2 of combined collecting area and with a spectral reso- 
lution given roughly by: 

AA s=a 0.12A x size of source in arcminutes/spectral order 

(2) 

for an extended source larger than 10 arcseconds. The wave- 
length band extends from 5 to 38 A, which samples Si K shell 
transitions to C K transitions. The field of view is effectively 
5 arcminutes by 1 degree, which is large enough to capture the 
entire cooling-flow region. The former is set by the width of 
the CCD array in the cross-dispersion direction, and the latter 
is set by self-vignetting of the telescope shells. For the cooling- 
flow clusters, data selection cuts are tailored to include mostly 
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emission from the cores of the clusters. The background is pro- 
duced primarily by a variable flux of soft protons spread nearly 
uniformly across the focal plane. There is no source-free re- 
gion for extended sources, so it has to be modeled. The cross- 
dispersion direction provides a one-dimensional image (11 arc- 
seconds FWHM) of the source, which we use to construct a 
model of the surface brightness of the source. 

4. SAMPLE SELECTION AND FLUXED SPECTRA 

Our sample of clusters (listed in Table 2) is a biased selec- 
tion of 14 compact clusters and groups of galaxies chosen to 
exploit the spectral sensitivity of the RGS. The sample includes 
clusters and groups at a range of temperature from 1 to 10 keV. 
All of these clusters and groups were expected to host cooling- 
flows ranging from 1 to 1000 solar masses per year (see Table 
3). 

To produce effective-area and exposure-corrected spectra of 
these objects, we adopt the following procedure: We process 
the events by a development version of the Science Analy- 
sis System (SAS) software (version 5.3), which accomplishes 
event reconstruction, aspect correction, CCD-pulseheight cor- 
rections, and focal plane reconstruction. We only use events 
in time intervals where the background was less than 6 counts 
per second in the 0.35 to 1.9 keV energy band. We then select 
photon events satisfying a 2 arcminute wide cross-dispersion 
cut. Events must also satisfy a first or second order joint 
dispersion-coordinate/CCD-pulseheight cut based on the char- 
acteristic resolution of the CCD, and additional broadening due 
to the 2 arcminute source extent. In this way, we select pho- 
tons roughly produced in a 2 arcminute square centered on each 
cluster. There is also a significant contribution of photons from 
outside this square, which is self-consistently in our modeling. 
Finally, a wavelength is assigned to each photon based on the 
nominal center of the cluster and the instrument bore-sight. 

In order to account for the instrument response and efficiency 
of these selection cuts, we use a Monte Carlo method, as dis- 
cussed in detail by Peterson, Jernigan & Kahn (2002). The 
Monte Carlo method is used to account for off-axis behav- 
ior of the response and the effect of arbitrary selection cuts 
and transformations made on the data. Additionally, it is re- 
quired in the astrophysical modeling as we discuss in §6. The 
Monte Carlo calculation contains all known effects of the mir- 
ror shells, grating arrays, and CCD response. It has calibra- 
tion limitations similar to the response in SAS version 5.3. 
The Monte Carlo method generates focal plane coordinates and 
CCD-pulseheights based on a model for the X-ray spatial and 
spectral distribution. 

We generate photons having a flat wavelength spectrum and 
surface brightness profile given by a modified (3 model (de- 
scribed in §6). The parameters are listed in Table 3. The Monte 
Carlo method uses focal plane maps of the exposure time pro- 
duced by the SAS to select detected events. The same selection 
cuts that are applied to the data, are applied to the simulated 
photons as well to account properly for the efficiency of the se- 
lection regions. Extracting the simulated photons as we did for 
the source photons produces a wavelength histogram in units of 
cm 2 s A. We then take the histogram of the source photons and 
subtract a set of background simulated events. The background 
model is described in detail in §6. The final histogram is di- 
vided by the exposure-area histogram and produces the fluxed 
spectra shown in Figures 3a, 3b, and 3c. The spectra are not 
corrected for absorption in the interstellar medium. 



5. QUALITATIVE RESULTS ON THE SAMPLE 

The spectra shown in the three panels of Figure 3 are sorted 
inversely by total luminosity, and therefore also roughly in- 
versely by temperature and M. Below, we discuss the spectra 
grouped by their outer temperature. 

7-10 keV Clusters: The RGS spectrum of Abell 1835 has 
already been presented in Peterson et al. (2001). The spec- 
trum exhibits intermediate temperature (kT e w 3 keV) plasma 
in addition to the 8 keV background plasma as shown by the 
detection of Fe XXIV and XXIII. O VIII Lya is clearly de- 
tected. Peterson et al. (2001) identified the gross inconsisten- 
cies between the measured spectrum of this source and the spec- 
trum predicted by the standard cooling-flow model. These con- 
straints are derived from the lack of Fe XVII-XXII. Abell 665 
has a similar spectrum, but the observation of that source was 
plagued by a very high background rate. 

4-7 keV Clusters: Abell 1795, Hydra A, Abell 496, and 
Abell 4059 are 4 to 7 keV clusters. Fe XXIV-XXII, O VIII, 
Mg XII, Ne X, and Si XIV emission lines are detected in all 
cases. The strong Fe XXIV-XXII complex indicates significant 
cooling down to kT e w 2 to 3 keV. There are no detections of Fe 
XVII-XXI, however, in clear contradiction to the predictions of 
the cooling-flow model for these systems. 

2-4 keV Clusters: The lowest temperature clusters, 
2A0335+096, Sersic 159-03, Abell 262, Abell 1837, Mkw3s, 
Abell 2052, and M87 have the bulk of their emission from the 
Fe L temperature range. These spectra show clear detections 
Fe XXIV-XIX, O VIII, Mg XII, Ne X, Si XIV. Various differ- 
ent ions produce stronger emission lines depending on the clus- 
ter temperature. Weak emission from the helium-like charge 
states are observed for silicon and magnesium, and in the best 
resolved spectra hydrogen-like C and N are detected. These 
spectra also generally have weak emission from Fe XVII and 
XVIII indicating some plasma below 1 keV, but less than would 
be predicted by the standard cooling-flow model. 

1 keV Groups: NGC 533 is the lowest temperature system 
and has emission from the same ions as in the 2-4 keV clusters. 
In this case, Fe XVII and XVIII are strong due to the lower 
temperature of the system. However, this spectrum is also in- 
consistent with the standard cooling-flow predictions, since Fe 
XVII is not much stronger than Fe XVIII and OVII is not de- 
tected. 

Composite Spectrum: As further confirmation of the real- 
ity of some of the detections discussed above, we generate a 
composite spectrum by combining all of the individual clus- 
ters. This composite is shown at the bottom of Figure 3. It was 
constructed by first shifting the wavelengths of photons by the 
cosmological redshift and then co-adding the counts-weighted 
spectra. The individual area-exposures were also added using 
the same technique. The counts-weighted spectrum was then 
divided by the sum of the area-exposures. We excluded M87 
from this analysis since it would otherwise dominate the result. 

In the composite spectrum, lines from Si XIV & XIII, Al 
XIII & XII, Mg XII & XI, all Fe L ions, Ne X, O VIII, C VI, 
and N VII are all clearly visible. No neutral O K edge is ob- 
served at the redshift of the clusters. This would be the most 
prominent absorption feature if a large amount of cold material 
with a large covering fraction was absorbing the flow. A feature 
which appears to be a blue-shifted O K absorption edge is due 
to galactic absorption, and its depth is consistent with the mean 
interstellar column density. 
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Several qualitative conclusions can be drawn from the obser- 
vations: First, in all cases, there is a deficit of emission expected 
from lower temperatures. Second, in most cases, a continuous 
distribution of temperatures is nevertheless required. Finally, 
absorption by intervening cool gas at the redshift of the cluster 
cannot explain the dearth of soft X-ray emission. In §6 below, 
we perform quantitive spectral fits to further explore the depar- 
tures from the cooling-flow model. 

6. FITTING METHODS AND MODEL 

In this section, we set quantitative limits on emission from 
various parts of the temperature distribution. To derive these 
limits, we use the Monte Carlo methods presented by Peterson, 
Jernigan & Kahn (2002), and implicitly used by Peterson et al. 
(2001) and Xu et al. (2002). The basic procedure is based on an 
astrophysical model for the spatial and spectral dependence of 
the emission. The Monte Carlo approach is required because a 
different spectrum at each projected spatial position is required 
to test the cooling-flow model in detail. We randomly generate 
photons with an associated dispersion angle, cross-dispersion 
angle, and CCD pulseheight value. The resulting count distri- 
butions are then compared with the raw data, after the various 
data selection cuts and transformations. In fitting for global 
parameters, such as the background normalization, we find the 
best fitting solution by using multivariate methods, as described 
in Peterson, Jernigan & Kahn (2002). For the abundances and 
temperatures, we use a \ 2 statistic of the combined extracted 
first and second order spectra for both instruments and we it- 
eratively adjust the surface brightness distribution to match the 
cross-dispersion profile. 

For the cluster emission we adopt a relatively simple model, 
so as to reduce sensitivity to fitting biases. For the surface 
brightness, we use a spherical (3 profile where the core radius is 
left free and the f3 parameter is fixed to the value determined 
from EPIC spectral fits (Kaastra et al. 2002). Additionally, 
we allow the normalization of the emission inside of a three- 
dimensional radius, r coo i, to be larger than the normalization 
outside of the radius. In this way, the spatial profile can be 
much more peaked than the standard (3 profile. This three di- 
mensional distribution is then projected on the sky. The precise 
shape of the spatial distribution is not critically important for 
determining the cooling luminosity limits as long as it repro- 
duces the rough behavior of the emission profile. 

Outside of r coo i, the emission is set to an isothermal temper- 
ature. Inside that radius, we fit for the normalization of the dif- 
ferential emission measure distribution below the upper temper- 
ature, Tq. The emission measure is divided into several temper- 
ature bins between jT and T , jTq and jT , |r and jTq, and 
-j^7b and |7b. This choice is arbitrary, but it provides robust fit- 
ting solutions since the fractional ionization curves are roughly 
equally spaced in the logarithm of the temperature. Preliminary 
fits indicate that only negligible hot ambient plasma coexists in- 
side the cooling radius, so additional isothermal emission inside 
this radius is ignored. Within each temperature bin, we use the 
isobaric radiative cooling-flow model to predict the shape of 
the emission measure distribution, but this is not critically im- 
portant since we do not have the spectral sensitivity to sample 
the emission measure distribution in very fine intervals. This 
approach clearly forces the coolest emission to conform to a 
specified spatial distribution. Further detailed analyses are re- 
quired to pinpoint its exact spatial location. Our spatial model, 
however, seems to be compatible qualitatively with the rela- 



tively small differences in the observed emission lines spatial 
profiles of M87 (Sakelliou et al. 2002). 

We assume uniform spatial abundances throughout the clus- 
ter emission. This considerably simplifies the fitting procedure, 
and makes the limits on the coolest emission conservative. It 
is known that abundance gradients exist in clusters, but they 
are generally flat within the cooling-flow region and weaken 
whenever more freedom is given to the temperature distribution 
(Molendi & Pizzolato 2001). The iron, neon, oxygen, magne- 
sium, and silicon abundances are left as free parameters. All 
other elements are tied to iron since they contribute few counts. 
In Abell 1835 and Abell 665, the magnesium, neon, and sili- 
con abundances are tied to iron since the emission is from very 
high temperatures. The absorption column density is left as a 
free parameter to account for variations along the line of sight 
to the cluster. In NGC 533 the absorption is set to the galactic 
value since it has little low energy continuum emission. We also 
ignore the effects of resonant scattering (cf. Xu et al. 2002), 
which redistribes emission line photons within our aperture, but 
otherwise results in the same detected emission line flux. 

For the background, we use a semi-empirical model cali- 
brated on blank sky Lockman Hole observations (XMM Rev- 
olutions 0070/0073). The model includes a spatial model for 
soft protons, low energy detector readout noise, and character- 
izations of the in-flight Al K and F K calibration sources. All 
parameters are frozen in this model except the relative normal- 
ization of the particle component, and the overall normalization 
which can vary by factors of 10 from observation to observa- 
tion. The background is relatively flat in wavelength. 

The model has the following free parameters: local column 
density, normalization of each part of the cooling flow emis- 
sion measure distribution, abundances of magnesium, neon, sil- 
icon, oxygen, and iron, background temperature (7b), particle 
background normalization, position of the source in the cross- 
dispersion direction, core radius, cooling radius, and overall 
normalization. For each spectrum, we apply selection cuts dis- 
cussed in §4. However, the background model parameters and 
position of the source are determined before applying any data 
selection cuts. 

In addition to the model described above, we also set lim- 
its on a model wherein an intrinsic absorber embedded in the 
cooling-flow volume is invoked to suppress the expected soft 
emission. If the absorber is evenly embedded, the transmission 
function is given by, (l -e" T(£) ) /t(E), where t(E) is the pho- 
toelectric optical depth as a function of energy, as opposed to 
the usual exponential form. For this case, we otherwise take 
the isobaric cooling-flow emission measure distribution, with- 
out allowing any variation from one temperature bin to the next. 

The errors are quoted at the 90% statistical confidence level 
for setting limits on the parameters. The uncertainty in the 
wavelength scale of 8 mA and the characteristic line spread 
function uncertainty of 5 mA make negligible contributions to 
these errors for extended sources like clusters. The effective 
area uncertainty is of order 10%. Additionally, there is known 
to be non-statistical noise at the 5% level of the source flux 
in the RGS spectra, caused by systematic dark current varia- 
tions which can sometimes produce false features. Either of 
these effects can produce correlated errors in adjacent wave- 
length channels. The global temperature, absorption column 
density, and overall normalization depend on many wavelength 
channels and are thus dominated by systematic uncertainties of 
this kind. These parameters are therefore never quoted with 



High Resolution X-ray Spectroscopic Constraints on Cooling-Flow Models 



5 



uncertainties below 200 eV, 5 x 10 19 cm~ 2 , or 10% in the nor- 
malization, respectively, which would all modify the spectrum 
characteristically by ~ 10%. 

We adopt the MEKAL plasma model (Mewe, Kaastra, & 
Liedahl 1995) as implemented in XSPEC vl 1 (Arnaud 1996) as 
the collisional ionization equilibrium model used to predict the 
spectrum. The absorption cross-sections are taken from Morri- 
son & McCammon (1983). Abundances are quoted relative to 
Anders & Grevesse (1989). Galactic absorption column densi- 
ties are compared to Dickey & Lockman (1990). We assume 
Ho = 70 km/s/Mpc, A = 0.7, O m = 0.3 throughout our analysis. 

Finally, all measured spectral mass deposition rates are com- 
pared to those derived in the corresponding EPIC data set. That 
analysis takes into account temperature and density gradients 
by solving the hydrodynamic equations, assuming only radia- 
tive cooling and mass drop out according to standard techniques 
(e.g. White, Jones, & Forman 1997). The mass deposition rates 
are interpolated at the radius for our r coo i parameter and are 
within uncertainties (typically 50%) of other published values. 
A more extensive discussion of the EPIC analysis will be pre- 
sented elsewhere (Kaastra et al. 2002). 

7. COOLING LUMINOSITY LIMITS 

Our derived spectral fits are shown in Figure 4a, 4b, and 4c 
and the cross-dispersion profiles are shown in Figure 5. The 
best fit parameters are given in Tables 3, 4 and 5. In each 
plot the red line shows the fitted empirical model. The green 
line shows the standard isobaric cooling-flow model without 
altering the temperature distribution and maintaining the same 
normalization. In almost all cases, the spectral features and 
the intrinsic line profiles are well-fit for the empirical model. 
Due to the low signal to noise of the observation of Abell 1837 
and Abell 665, we could not place strong constraints on the 
cooling-flow model for those clusters. The other clusters, how- 
ever, show clear deviations from the cooling-flow model. 

There are two significant residuals that should be noted. The 
long wavelength region of the spectrum is not perfectly fit. We 
believe this is due to errors in the effective area calibration, im- 
perfect background modeling, and spatially broadened (~ 3 A) 
emission lines of O VIII, O VII, N VII, and C VI emission 
lines from the diffuse soft X-ray background (see McCammon 
et al. 2002). These discrepancies have a minimal effect on our 
conclusions about the temperature distribution, however, since 
those rely primarily on the Fe L region. Another problem is 
evident in the detailed line fitting of the spectrum of NGC 533. 
This is probably due to an underprediction of the 2p-3s lines of 
Fe XVII and Fe XVIII (at 16 and 17 A, respectively). There is 
a known problem with the excitation rates for these transitions 
(Beiersdorfer et al. 2002), as has been discussed for the case of 
NGC 4636 by Xu et al. (2002) and Capella by Behar, Cottam, 
& Kahn (2001). 

The differential luminosity distribution for all clusters is plot- 
ted in Figure 6. The luminosity is normalized with respect to 
the prediction from the isobaric radiative multi -phase model us- 
ing Equation (1) and M from Table 4. The presence of points 
above the dotted line is not unexpected, since these can result 
from the background uncooled hot plasma. This can also re- 
sult from the ambiguity in defining a distinct cooling radius. 
However, the data should closely track the line y = 1 at lower 
temperatures, if the cooling-flow model is correct. Many upper 
limits, however, are well below that line at temperatures near 
a third of the background temperatures. As can be seen, the 



well-measured clusters show clear deviations from isothermal- 
ity with significant plasma existing down to one-quarter of the 
background temperature. In Figure 7, we plot the same results 
with the horizontal axis now scaled to the background temper- 
ature. Here the discrepancies from the cooling-flow prediction 
look more systematic. In particular, the differential luminosity 
distribution appears to be roughly consistent with the expres- 
sion, 



dL 5 Mk ( T 

al 2 [im p \ 1q 



(3) 



where a w 1 to 2 instead of 0, as expected for the isobaric ra- 
diative cooling-flow model. There is still significant scatter in 
the normalization and slope of this relation, which could be due 
to various systematic effects, or to real differences between the 
cooling-flows. Further detailed analyses that start with the as- 
sumption of Equation (3) might yield further insight into the 
scatter. However, it is clear that no cluster has a temperature 
distribution consistent with the a = case, and that in each case 
the failure of the cooling-flow model occurs at a fraction of the 
virial temperature rather than at a fixed value. 

The derived abundances as a function of ambient temperature 
are plotted in Figure 8. The abundance of iron declines slightly 
with more massive clusters as indicated in earlier ASCA obser- 
vations (Mushotzky et al. 1996, Fukazawaet al. 1998), although 
the abundance measured here is somewhat higher. The emis- 
sion weighted iron abundance ranges from one third to 70% 
with higher values in lower mass clusters. This presumably 
reflects the increase in stellar mass fraction for lower mass sys- 
tems. There are no obvious trends in the other abundance ratios 
and there is considerable scatter in the points. This is possibly 
caused by subtle fitting biases. In particular, the abundances 
reflect the emission-weighted abundance average in an aper- 
ature that varies from cluster to cluster and there is often an 
anti-correlation between the absolute abundances and the width 
of the differential emission measure distribution. Additionally, 
abundances of elements other than iron are measured off the 
peak of their fractional ionic abundance, and depend sensitively 
on the temperature distribution. 

Nevertheless, all of the abundances seem to cluster around 
particular ratios. The oxygen to iron ratio is 70 ± 20 % (1 a) 
of the solar value. The magnesium to iron ratio is 100 ± 40 %, 
the neon to iron ratio is 110 ± 40 %, silicon to iron is 230 ± 
80%. Abundance patterns for at least these four elements re- 
semble those found in Xu et al. (2002) suggesting a common 
enrichment history between ellipticals and clusters. These par- 
ticular ratios do not appear to fit into a simple superposition of 
Type la yields plus Type II supernovae integrated over a single 
initial mass function, as in Gibson et al. (1997). They more 
closely resemble low-metallicity high-mass supernovae yields, 
but further analyses with more attention to additional elements 
are needed to establish the pattern completely and determine 
the spatial distribution. Additionally, carbon and nitrogen abun- 
dances could be roughly estimated in M87 and Abell 496 to be 
C/O ~ 2 and N/O ~ 2.5 with large uncertainties. This implies, 
however, a non-negligible fraction of chemical enrichment has 
been processed through stellar winds. Further detailed analyses 
of the measured abundances and their spatial distributions will 
be presented elsewhere (Tamura et al. 2003). 

The derived intervening column densities are generally con- 
sistent with the expected galactic values. The only large de- 
viation is with 2A0335+096, which has a high galactic col- 
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umn density and the measurements could be more affected by 
systematic uncertainties. In addition, allowing an intrinsic ab- 
sorber with the standard cooling-flow model produced worse 
fits than the standard cooling-flow model alone. Although the 
absorber can remove significant cold emission with column 
densities in excess of 2 x 10 21 cm" 2 , it also removes the low 
energy tail of the bremsstrahlung due to the core, which pro- 
vides a significant fraction of the continuum. Limits on in- 
trisic absorbers are given in Table 4 and are generally less than 
2 x 10 20 cm- 2 . 

8. DISCUSSION 

In §7, we demonstrated the presence of significant quanti- 
ties of plasma just below the ambient temperature, 7b, down to 
To /2 and the clear failure of the cooling-flow model at lower 
temperatures. We also showed the consistency of the sample 
with the empirical differential luminosity distribution given by 
Equation (3) with the parameter, a, equal to » 1 to 2. Here, 
we consider several ideas which have been suggested to ex- 
plain these results. We first discuss mechanisms which mod- 
ify the isobaric cooling-flow model, but involve no real mod- 
ification to the general paradigm of radiatively-driven flows in 
gravitationally-relaxed cluster cores. Then, we discuss substan- 
tial modifications to the general cooling-flow process by inclu- 
sion of additional heating mechanisms or non X-ray cooling 
channels. 

8.1. Modifications to the Standard Spectral Prescription 

We have compared the data to the standard isobaric cooling- 
flow model. A modification to this model is expected if a weak 
magnetic field is amplified as cooling plasma compresses in or- 
der to maintain pressure equilibrium. At the lowest temper- 
atures then, these blobs cool isochorically. This would only 
modify the expression in Equation (1) by changing the | to |, 
which is still inconsistent with the data displayed in Figures 6 
and 7. 

The compression of the gas due to the dark matter-dominated 
gravitational potential is included in our calculation of M by 
measuring directly both the temperature and density spatial gra- 
dients, and has generally been included in previous work (e.g. 
Allen 2000, Peres et al. 1998). This is important if much of the 
gas is actually flowing, and the effect on the temperature dis- 
tribution is given in Nulsen (1998). It has a relatively small 
effect on the spectra apart from the normalization, however, 
and should not affect the luminosity at the lowest temperatures 
since most of the plasma is thought to drop out of the ambient 
medium before it flows (Johnstone et al. 1992). Even if there is 
no mass drop out, however, the | in Equation (1), would only 
be replaced by | - X T / X p where X T is the radial logarithmic 
temperature gradient and X p is the radial logarithmic density 
gradient. X T / X p is always observed to be negative in the cores 
of clusters. 

Two recent explanations to the soft X-ray cooling-flow prob- 
lem have been proposed that involve no real modification to 
the multi-phase model in Equation (1). First, Peterson et al. 
(2001) suggested that differential absorption by colder mate- 
rial could absorb more of the emission from cold ions. This 
is somewhat implausible, since one would expect hot and cold 
parcels of plasma to have roughly equal probability of being 
obscured. Nevertheless, the result could be interpreted in the 
context of this model if the effective absorption column density 
scaled as a strong function of the occupied volume. This would 



require the absorber to have an almost identical spatial distri- 
bution to the currently cooling parcels of gas from the cooling- 
flow and a small filling factor. It is impossible to rule out this 
model in its most extreme form, since an infinite column den- 
sity selectively applied to emission at a given temperature is, of 
course, identical to no emission at that temperature. The most 
model-independent limits against this explanation comes from 
Bohringer et al. (2002) looking at absorption in spectra of AGN 
in clusters. 

An explanation offered by Fabian et al. (2001) and developed 
in Morris & Fabian (2002) involved a locally non-uniform dis- 
tribution of metals in the ICM, which modifies the spectrum 
significantly. This requires that the metal abundance variations 
occur on scales smaller than cooling parcels of plasma, and that 
the variations sample a particular bimodal distribution. In its 
simplest form, this model would predict that the cooling-flow 
predictions would exceed the measured spectra near roughly 
a common temperature near 1 keV, where line cooling domi- 
nates over bremsstrahlung for super-solar abundances. That is 
not what is observed, however, and this model is therefore not 
compatible with our results for the low temperature groups. 

8.2. Modifications to Radiatively -Driven Flows 

A variety of additional physical processes have been pro- 
posed to operate in the cores of clusters, which could, in princi- 
ple, account for the soft X-ray cooling-flow problem (Peterson 
et al. 2001, Fabian et al. 2001, Bohringer et al. 2002, Sasaki 
& Yamasaki 2002, and Brighenti & Mathews 2002 and refer- 
ences therein). These include time-dependent heating by AGN 
outflows (e.g. Rosner & Tucker 1989, Tabor & Binney 1993, 
David et al. 2001), electron thermal conduction from the outer 
regions of clusters (e.g. Tucker & Rosner 1983, Stewart et al. 
1984, Bertschinger & Meiksin 1986, Voigt et al. 2002, Zakam- 
ska & Narayan 2001, Ruszkowski & Begelman 2002), contin- 
ual sub-cluster mergers (Markevitch 2001), and interactions of 
the ICM with dark matter (Qin & Wu 2001). Others invoke 
rapid cooling induced by mixing at interfaces with cold clouds 
(Begelman & Fabian 1990, Fabian et al. 2001), dust mixing 
(Fabian et al. 2001), or acceleration of relativistic particles. 

Most calculations have succeeded in establishing that the 
above mechanisms are energetically important, i.e these mech- 
anisms have total heating powers that can roughly cancel ra- 
diative cooling energy losses, or that alternative cooling chan- 
nels have sufficient luminosity to cool the plasma through other 
means. The total energy balance (i.e. the sum of cooling lu- 
minosity minus heating power) determines the total mass de- 
posited from the cooling-flow, and thus directly addresses the 
classic cooling-flow problem. The total energy balance, how- 
ever, only affects the soft X-ray cooling-flow problem by re- 
ducing the total normalization of Equation (3) by a fixed value. 
It does not change the predicted differential luminosity distri- 
bution. 

The critical difference, however, between the soft X-ray 
cooling-flow problem and the classic cooling-flow problem is 
that the latter requires a clear explanation for why X-ray cool- 
ing does not appear to be carried to completion. It is difficult to 
find a relevant heating, mixing, or cooling time scale that would 
be comparable to the X-ray cooling time. If the time scale for 
a given process is too short it will overwhelm the cooling-flow 
(a = oo), and if the time-scale is too long, it is dynamically 
unimportant (a = 0). Below, we discuss the proposed physical 
processes and whether they can account, both energetically and 
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dynamically, for the missing soft X-ray luminosity. 

8.2. 1 . External Heat Sources 

There are three requirements for additional heating mecha- 
nisms to be compatible with the observations. First, the total 
time averaged heating power, <P>, has to roughly cancel the 
radiative losses so that the expression, 

<P>^p (4) 

is near unity. The denominator in the expression varies among 
the clusters in our sample by four orders of magnitude, so the 
proposed heating process must operate on a variety of scales. 
Second, the heating has to be distributed spatially throughout 
the cooling-flow volume to cancel cooling everywhere. Third, 
the process has to be self-regulating, so that the time scale for 
heating remains comparable to the cooling time for all clusters. 

AGN outflows: Time-dependent AGN outflow heating mod- 
els have been considered by a variety of authors (e.g. Ros- 
ner & Tucker 1989, Tabor & Binney 1993, Churazov et al. 
2001, Bruggen & Kaiser 2001, Quilis 2001, David et al. 2001, 
Nulsen 2002). Buoyant bubbles carrying relativistic plasma ap- 
pear to be a common phenomena in clusters with central AGNs. 
The thermal energy seems to be enough to heat cooling-flows 
through cosmic ray interactions and mechanical heating, but it 
is unclear whether this energy gets properly channeled into the 
cooling volume (e.g. Loewenstein, Zweibel, Begelman 1991, 
Fabian et al. 2001). Note that it is essential that the heat be 
distributed evenly throughout the region which is thermally un- 
stable. In addition these models must be made self -regulating to 
counteract cooling at a rate proportional to the mass deposition 
rate, and with periods of heating roughly as long as periods of 
cooling. Clearly, it requires a significant degree of fine-tuning. 

Electron Thermal Conduction to the Background 
Plasma: There is considerable thermal energy in the outer 
regions of clusters that can destroy any existing cooling-flow 
through electron thermal conduction (e.g. Tucker & Rosner 
1983, Stewart et al. 1984, Bertschinger & Meiksin 1986). The 
size of cooling-flows are only a few electron mean free paths 
in the absence of magnetic fields. The critical question is to 
what level is conduction suppressed by tangled magnetic fields, 
an issue which continues to be debated theoretically (Chandran 
& Cowley 1998, Narayan & Medvedev 2001). Observation- 
ally, conduction is suppressed by factors near 100 in identi- 
fied cold fronts (Ettori et al. 2002, Markevitch 2000, Vikhlinin, 
Markevitch & Murray 2001). Voigt et al. (2002), Zakamska & 
Narayan (2001), Fabian, Voigt, & Morris (2002) have demon- 
strated that the heat flow from the outer regions of clusters with 
a small suppression (> 10" 1 ) in the Spitzer conductivity would 
appear to cancel radiative losses in many clusters. The spatial 
distribution of the heating and overall energetic requirements 
appear to be satisfied by conduction models, but there is no 
explanation for why the cluster would cool to their current tem- 
perature distribution. Since conduction suppresses temperature 
gradients by definition, this mechanism alone does not solve the 
dynamical problem presented here. 

Sub-cluster Mergers: Markevitch (2001) has shown that 
clusters previously thought to be fully relaxed exhibit tempera- 
ture fronts consistent with plasma exhibiting large bulk motions 
in the gravitational potential. This provides another source of 
energy that has not been dissipated and therefore leads to an 
increase of the cooling time above previous estimates (Gomez 



et al. 2002), so that the time scales for radiative cooling and 
dynamical relaxation may be similar. There is no reason for 
this process to be self-regulating, however, and we would ex- 
pect a much larger difference in the temperature distributions, 
depending on each cluster's merger history. This explanation 
also requires a conspiracy of factors to both cancel radiative 
cooling in global energetics, and ensure that the mergers occur 
with a frequency that allows some cooling. Future numerical 
simulations may test this further, and presumably observations 
of cooling-flows at a different epoch would not show the same 
effects. 

8.2.2. Rapid Cooling Mechanisms (Heat Sinks) 

The radiative isobaric cooling-flow model assumes that all 
of the thermal energy is released in X-rays at high tempera- 
tures. There may, however, be additional contributions from 
other cooling processes. There are three main requirements for 
additional cooling channels to explain the observations, which 
are similar to, but slightly different from the heating require- 
ments. The first is that the total power in the coolant be com- 
parable to the missing soft X-ray luminosity. Any coolant with 
power, Pcooiam, reduces the total X-ray emission by a factor, 



1 _i_ 'cool (lilt 

The second requirement is that the ratio P ' Y £'""' needs to have 
a temperature dependence consistent with Equation (3) or that 
the cooling channel is self-regulating in the same sense as the 
discussed heating models. The third requirement is that the en- 
ergy should be released with a similar spatial distribution to the 
lowest temperature X-rays. 

UV cloud interfaces: Begelman & Fabian (1990) and 
Fabian et al. (2001) discussed the possibility that hot electrons 
are cooled conductively by interfaces with cold clouds. This 
leads to emission in the UV where the cooling function is the 
highest, and is consistent energetically with the large observed 
Ha luminosities (Heckman et al. 1989, Crawford et al. 1999) 
of 10 42 to 10 44 ergs/s in cooling-flows. However, by itself, this 
model does not explain the observed temperature distribution in 
the soft X-ray band, since high temperature electrons are cooled 
by this process as well. In addition, in such a picture highly 
charged ions should also impact the cloud interfaces, result- 
ing in charge exchange, which produces copious soft X-ray line 
emission. The spatial distribution of Ha is remarkably similar 
to the coolest X-rays (Ettori et al. 2002) and the total luminosity 
is marginally sufficient to account for the missing soft X-ray lu- 
minosity (Fabian et al. 2001), but again the dynamical problem 
of partial cooling is not easily solved by this mechanism alone. 

Dust Mixing: Some clusters have been established as strong 
IR sources, typically of order 10 43 to 10 45 ergs/s (Edge et al. 
1999). Fabian et al. (2002) demonstrated that in all systems 
studied, the missing thermal energy in soft X-rays is more than 
accounted for by the IR luminosity. Hot electrons impinging 
on dust grains could collisionally heat the grains and this in 
turn cools the hot electrons. The cooling function for this pro- 
cess has been discussed in Burke & Silk (1973) and Draine & 
Salpeter (1979) and is 100 times larger than the X-ray cooling 
function for a typical ISM composition of dust grains. As such, 
this explanation appears promising, but there is no strong tem- 
perature dependence to either the X-ray or dust cooling function 
that would systematically reduce the overall level of soft X-ray 
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emission. Presumably, the success of this process depends on 
magnetic field topologies and a steady-state process of destruc- 
tion and creation of dust grains. Suppression of conduction by 
magnetic structures in some regions could play a critical role in 
regulating both UV cooling and dust mixing, and fast magnetic 
reconnection could in turn realign those structures. The dynam- 
ical problem remains unsolved here too, and if dust cooling is 
important for cluster X-ray plasmas, it creates an even larger 
cooling-flow problem. 

Relativistic Particles and Alfven Waves: The increase in 
magnetic field density as cooling blobs compress may allow for 
considerable energy to be released via Alfven Waves or mildly 
relativistic particles. In fact, the steeper synchrotron indices 
and short particle lifetimes in cluster radio mini-haloes have 
been invoked to argue for reacceleration of particles in cooling- 
flows (e.g. Petrosian 2001, Gitti, Brunetti, & Setti 2002). This 
process should, however, only liberate a small fraction of the 
thermal energy and it provides no natural explanation for the 
observed temperature distribution. 

8.3. General Considerations and Future Work 

Implicit in our discussion above is that only one mechanism 
explains the observations. Given that many of the above ideas 
are energetically viable, that may be unlikely. It may be that 
the soft X-ray cooling-flow problem and the classic cooling- 
flow problem are solved by different solutions. It may also be 
possible that a combination of processes could have different 



effects. For example, heating cooling-flows both from the in- 
side with an AGN and from the outside through conduction, 
results in a more stable heating method than each mechanism 
does individually (Ruszkowski & Begelman 2002). 

We have demonstrated that the soft X-ray cooling-flow prob- 
lem is not an isolated phenomenon and appears on a range of 
scales. The empirical characterization offered here may provide 
clues to the ultimate solution. Clearly, the model for the mul- 
tiphase distribution is inadequate. More detailed comparisons 
between observations at X-ray and other wavelengths may more 
carefully address which of the physical mechanisms are impor- 
tant energetically. Much more theoretical work is needed in ad- 
dressing why any of the proposed mechanisms would dynami- 
cally conspire to produce the observed temperature distribution. 
These descriptions, however, are just beginning to be developed 
and it remains a theoretical challenge either to cool keV plasma 
with little resulting soft X-ray radiation, or to balance radiative 
cooling both energetically and dynamically. 

This work is based on observations obtained with XMM- 
Newton, an ESA science mission with instruments and con- 
tributions directly funded by ESA Member States and the US 
(NASA). Work on the RGS at Columbia University and U. C. 
Berkeley is supported by NASA. The laboratory for Space Re- 
search, Utrecht is supported by NWO, the Netherlands Organi- 
zation for Scientific Research. JRP acknowledges many helpful 
conversations with other members of the Columbia team. 
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Important line blends in Cluster Soft X-ray Spectra 



Ion 


Wavelengths 


Temperature 


Ion 


Wavelengths 


Temperature 




A 


keV 










Fe XXIV 


10.6, 11.2 


0.9 -> 4.0 


O VIII 


19.0, 16.0 


>0.2 




Fe XXIII 


11.0, 11.4 


0.8 -> 2.0 


O VII 


21.6, 22.0 


0.1 -> 


0.2 




12.2 




Si XIV 


6.2 


>1.0 




FeXXII 


11.8, 12.2 


0.6 -> 1.5 


Si XIII 


6.6,6.7 


0.2^ 


1.0 


FeXXI 


12.2, 12.8 


0.5 -> 1.0 


Al XIII 


7.2 


>1.2 




FeXX 


12.8, 13.5 


0.4^ 1.0 


Al XII 


7.8,7.9 


0.3 -> 


1.2 


FeXIX 


13.5, 12.8 


0.3 -> 0.9 


Mg XII 


8.4 


>0.7 




Fe XVIII 


14.2, 16.0 


0.3 -> 0.8 


Mg XI 


9.2,9.3 


0.1 -> 


0.6 


Fe XVII 


15.0, 17.1 


0.2 -> 0.6 


NeX 


12.2 


>0.4 






15.3, 16.8 




NelX 


13.5, 13.7 


0.1 -» 


0.3 


N VII 


24.8 


>0.1 










C VI 


33.7 


>0.1 











Table 2 

Basic Properties of the Sample 



Cluster 


XMM Rev. 


Eff. Exposure 


Redshift 


Lum. Distance 


Ang. Scale 






ks 




Mpc 


kpc/arcsec 


A 1835 


0101 


36 


0.2528 


1598 


6.2 


A 665 


0242 


20 


0.1818 


1042 


4.3 


A 1795 


0100 


40 


0.0622 


296 


1.4 


Hydra A 


0183 


38 


0.0538 


253 


1.2 


Ser 159-03 


0077 


38 


0.0580 


274 


1.3 


2A0335+096 


0215 


26 


0.0347 


158 


0.7 


A 4059 


0176 


54 


0.0460 


213 


1.0 


A 496 


0211 


29 


0.0328 


149 


0.7 


MKW3s 


0129 


39 


0.0442 


204 


1.0 


A 2052 


0128 


33 


0.0353 


161 


0.8 


A 262 


0203 


36 


0.0163 


71 


0.3 


A 1837 


0200 


50 


0.0372 


170 


0.8 


M87 


0097 


42 


0.0043 


19 


0.1 


NGC 533 


0195 


48 


0.0185 


82 


0.4 



Table 3 



Spatial Parameters 



Cluster 





^core 


fcool 






arcsec 


arcsec 


A 1835 


0.86 


79 


25 


A 665 


0.67 


53 


55 


A 1795 


0.74 


147 


55 


Hydra A 


0.73 


138 


44 


Ser 159-03 


0.79 


72 


29 


2A0335+096 


0.79 


156 


63 


A 4059 


0.78 


72 


156 


A 496 


0.64 


380 


72 


MKW3s 


0.72 


111 


25 


A 2052 


0.67 


258 


51 


A 262 


0.51 


163 


31 


A 1837 


0.67 


134 


12 


M87 


0.51 


178 


35 


NGC 533 


0.58 


119 


18 
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Table 4 

Abundances and Absorption Column Densities 



Cluster 


N H 


-^Galactic 


-^Intrinsic 





Ne 


Mg 


Si 


Fe 




10 20 cm- 2 


10 20 cm- 2 


10 20 cm- 2 













A 1835 


3±0.5 


2.32 


<6 


0. 


.19±0. 


.08 






















0. 


27 


A 665 


9±3.0 


4.24 






























0. 


33 


A 1795 


1±0.5 


1.19 


<2 


0. 


,38±0 


.11 


0. 


80±0. 


20 


0. 


,75±0. 


.25 


1 


,8±0. 


,6 





,42±0 


.04 


Hydra A 


4±0.5 


4.94 


<2 


0. 


,34±0. 


.06 


0. 


36±0. 


12 


0. 


,26±0. 


.13 


1 


,0±0 


.5 





,49±0. 


,08 


Ser 159-03 


0.5±0.5 


1.79 


<2 


0. 


,20±0. 


.07 


0. 


54±0. 


09 


0. 


,41 ±0, 


,19 


1 


,4±0. 


,3 





,55±0. 


.05 


2A0335+096 


28±0.5 


17.6 


<25 


0. 


53±0. 


.06 


0. 


62±0. 


10 


0. 


,50±0. 


,18 





,9±0. 


.2 





,56±0 


,05 


A 4059 


2±0.5 


1.10 


<4 


0. 


,52±0. 


,14 


0. 


70±0. 


18 


0. 


98±0. 


,22 


1 


,9±0 


.5 





,66±0 


,07 


A 496 


8±0.5 


4.58 


<5 


0. 


,44±0. 


,07 


0. 


70±0. 


22 


0. 


,68±0. 


,25 


1 


,9±0. 


.5 





,78±0. 


,08 


MKW3s 


2.5±0.5 


3.03 


<2 


0. 


,35±0. 


.06 


0. 


60±0. 


16 


0. 


,63±0. 


,16 





,6±0. 


.4 





,58±0. 


.06 


A 2052 


3±0.5 


2.73 


<3 


0. 


,46±0. 


,10 


0. 


,51±0. 


15 


0. 


84±0. 


,28 


1 


,4±0. 


.4 





,61±0, 


.06 


A 262 


9.5±0.5 


5.37 


<12 


0. 


.51 ±0 


.09 


0. 


37±0. 


22 


0. 


87±0. 


,23 


1 


,8±0. 


.4 





,68±0, 


,07 


A 1837 


4±0.5 


4.30 


































M87 


4±0.5 


2.54 


<2 


0. 


,40±0 


.04 


0. 


70±0. 


07 


0. 


,55±0. 


,06 





,9±0. 


1 





,54±0, 


,05 


NGC 533 


3.10 


3.10 


<18 


0. 


,41 ±0. 


.05 


0. 


,52±0. 


.25 


0. 


,98±0. 


.30 


2 


,0±0. 


.5 





,89±0. 


,15 



Table 5 
Cooling Luminosity 



Cluster 


16 8 


Mi . 

8 4 


Mi i 

4 2 


M, , 

2 1 


t isothermal 


kT 


^morphological 




Moyr 1 


Moyr" 1 


Moyr 1 


M Q yr 1 


ergs s" 1 


keV 


Moyr 1 


A 1835 


<200 


<300 


800±200 


5800±800 


6xl0 4S 


9.5±0.5 


1000 


A 665 










2xl0 45 


7.7±0.5 




A 1795 


<30 


<30 


<80 


380±40 


2xl0 45 


5.5±0.5 


300 


Hydra A 


<100 


35±20 


120±60 


180±50 


7xl0 44 


6.0±0.3 


180 


Ser 159-03 


<30 


<30 


<60 


210±30 


3xl0 44 


3.8±0.3 


79 


2A0335+096 


<84 


20±10 


40±20 


420±50 


4xl0 44 


3.2±0.3 


170 


A 4059 


<10 


10±5 


40±20 


100±10 


3xl0 44 


6.0±0.3 


60 


A 496 


<10 


<15 


25±10 


120±20 


8xl0 44 


4.7±0.3 


72 


MKW3s 


<64 


<10 


<20 


45±10 


3xl0 44 


3.7±0.3 


13 


A 2052 


<10 


<10 


15±5 


100±20 


4xl0 44 


3.4±0.3 


50 


A 262 


<10 


<2 


1.6±0.5 


10±1 


9xl0 43 


2.1±0.2 


2.0 


A 1837 










8xl0 43 


4.2±0.4 




M87 


<12 


<0.6 


0.6±0.2 


5.9±0.5 


4xl0 43 


2.0±0.1 


2.4 


NGC 533 


<5 


<1.6 


2.3±0.4 


5±0.5 


6xl0 42 


1.5±0.1 


1.5 
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Isobaric Multiphase Cooling-Flow Model 



a 
o 

o 
A 

o 

10 

O 




20 25 
Wavelength (A) 



FIG. 1. — The predicted spectrum of the isobaric multiphase model for a maximum temperature of 6 keV. Note that Fe L lines and O VIII are particularly 
prominent. 



0.250 



Isobaric Multiphase Cooling-Flow Model 



| 0.125 

S3 
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20 25 
Wavelength (A) 



FIG. 2. — The same plot as above but divided into different temperature ranges: 6 to 3 keV (red), 3 to 1.5 keV (yellow), 1.5 to 0.75 keV (green), and 0.75 to 0.375 
keV (blue). Note that the Fe L ions provide critical diagnostics of the temperature distribution; whereas, the Lyo transition from hydrogenic ions is produced at a 
large range in temperatures. 
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20 25 
Rest Wavelength (A) 



FIG. 3. — (Three panels) Fluxed and background subtracted spectra of the sample of clusters. Both RGS instruments and both spectral orders are included in all 
the plots. Prominent detected or expected emission lines are labeled (also see Table 1) as well as the neutral O K edge. All spectra are deredshifted so the horizontal 
axis is the wavelength in the cluster rest frame. The last panel shows all the cluster spectra added together with the exception of M87. The spectra are divided by 
the effective area so some regions of the spectrum contain more statistical noise than others. In particular, the long wavelength region and the regions between 20 
and 24 A and 10.5 and 14 A (in the lab frame) have lower effective area. 



Peterson et al. 




10 15 20 25 30 35 

Rest Wavelength (A) 



18 



Peterson et al. 




10 15 20 25 30 35 

Rest Wavelength (A) 



FIG. 4. — (Three panels) Comparison of the data (blue) with 1 a error bars, the empirical best fit model (red), and the standard cooling-flow model (green). 
Both the model and the data are fluxed and background subtracted. The cooling-flow model is not a best fit, but calculated by merely taking the soft X-ray flux in 
the empirical model and using the standard isobaric temperature distribution. Evident in all spectra is the severe overprediction of emission lines from the lowest 
temperatures. 
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Abell 2052 
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FIG. 5. — Cross-dispersion profiles (Id images) of each individual cluster and all profiles added together in the composite spectrum. The data are shown as a blue 
histogram and the model is the red line. Most discrepancies seem to result from the assumption of spherical symmetry in the model. 
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FIG. 6. — Differential Luminosity vs. Temperature computed by dividing the upper limit on M by the estimate derived from the imaging data. Each color 
represents a different cluster and lines connect the points, which are not upper limits. Points above the line are consistent with ambient plasma which has not cooled, 
but the isobaric multiphase model would predict that emission should closely track the dashed line at y = 1 or the line y = | for isochoric cooling. A large number 
of upper limits and points lie below those lines, however. The spectra generally become inconsistent with the model near | of the ambient temperature. 
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FIG. 7. — The same as the previous figure except plotted as function of the fraction of the ambient temperature. A general trend in the points which are not upper 
limits can be seen which follows a temperature distribution which is proportional to the fractional temperature or the fractional temperature to the second power 
(dotted lines). 
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FIG. 8. — Abundances and abundance ratios plotted at a function of the ambient cluster temperature. The iron abundance tends to increase for lower mass systems. 
No other obvious trends can be seen. The dashed lines are the average value for the sample and the 1 a ranges. 



